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Abstract 

We analyze the dynamics of concentrated polymer solutions modeled by a 2D Smolu- 
chowski equation. We describe the long time behavior of the polymer suspensions in a 
fluid. 

When the flow influence is neglected the equation has a gradient structure. The 
presence of a simple flow introduces significant structural changes in the dynamics. We 
study the case of an externally imposed fiow with homogeneous gradient. We show that 
the equation is still dissipative but new phenomena appear. The dynamics depend on 
both the concentration intensity and the structure of the flow. In certain limit cases the 
equation has a gradient structure, in an appropriate reference frame, and the solutions 
evolve to either a steady state or a tumbling wave. For small perturbations of the 
gradient structure we show that some features of the gradient dynamics survive: for 
small concentrations the solutions evolve in the long time limit to a steady state and 
for high concentrations there is a tumbling wave. 

AMS subject classification 35Q, 76A05 

1 Introduction 

In this paper we study qualitative properties of a Smoluchowski equation describing the 
dynamics of non-Newtonian complex fluids containing liquid crystalline polymers, in a 
concentrated regime. 

The model we use was introduced by M. Doi in ([7j) (see also jH|)- It identifles the 
polymers with inflexible rods whose thickness is much smaller than their length. We 
study here the case in which the fluid is two dimensional. This represents a simpliflca- 
tion that preserves many of the qualitative features of the physical three dimensional 
phenomenon (|22|)- This model as well as its three dimensional analogue has attracted 
much interest in the recent years 0, El , EH , [HI , [IS] , d , [H , Hi , EH , [Ml ) • 

The local probability measure associated with the polymers is of the form f{t, x, 9)d9. 
Here t is the time coordinate, x E denotes the spatial coordinate and 9 G [0, 27r] is a 
direction on the unit circle. The measure fd9 represents the time-dependent probability 
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that a rod with center of mass at x has an axis of direction 9 in the area element dO. 
The equation we study, in non-dimensional form, is: 



dtf + u- V,./ + de{Vf + deilCf) ■ f) = deef 



(1) 



where u is the velocity of the underlying fluid. For V^ii 



(u;)ij=i,2, 



V{x,9) = —u\ (x) cos 6 sin 6 — u\ (x) sin^ 6 + u\ (x) cos^ 6 + (x) sin 6 cos 6 (2) 

denotes the projection of V xU ■ (cos(0), sin(0))* on the tangent space in (cos(0), sin(0))*. 
This term describes the way the fluid influences the evolution of /. 

The term /C/ - the excluded volume potential, which accounts for the interaction 
between different rods - is given by 



where /c is a smooth function and 6 is a non-dimensional parameter measuring the 
concentration of the polymers in the fluid. Moreover 



ICf is the so called Maier-Saupe potential. This potential has been frequently used in 



The fluid velocity u obeys the Stokes or the Navier-Stokes equations forced by an 
appropriate average of / (|HI)- 

The rich dynamical behavior of the system poses significant numerical and analytical 
challenges (123; |2E])' We consider two levels of complexity: 

On a first level one neglects the influence of the fluid. This situation has been 
analyzed in [E])- It was shown that the system has a gradient structure and it is 
dissipative, i.e. the solutions starting from an arbitrary initial data end up in a fixed 
ball. We refine this analysis by showing that the a;-limit set of any solution consists of 
steady states. We also show that cj-limit set is reduced to only one steady state if some 
additional symmetry constraints are imposed. 

In the presence of the flow the dynamics become very complex, even when the 
flow has a simple structure. We consider the case of an externally imposed fiow with 
homogeneous gradient, i.e. the matrix is a given constant matrix. This is a situation 
of physical importance because it captures the local behavior of steady, smooth flows, 
and hence it is encountered frequently in the rheological literature (see P, chapter 5). 
A similar situation, when the fluid is a shear flow and has small amplitude, has been 
recently analyzed in [33] . 

We prove that the equation is dissipative in this case as well. It turns out that 
the dynamics depend strongly on both the intensity of the concentration b and the flow 
structure. Let us observe that in the case of homogeneous gradient flow V is independent 
of X and we can write 

V = uj scos(26' + a) 

with w G M, s > uniquely determined and a determined modulo 2tt. One can easily 
check that u) = "^2"^ so is the vorticity of the imposed flow. We can assume without 




(3) 




the literature (0, 0, 0)- 
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loss of generality that a = (if a / then f{t,6) = f{t,6 — a/2) satisfies the same 
equation but with cos(20 + a) m. V replaced by cos(26')). The parameters uj and s 
determine the flow structure. 

In certain limit cases the equation preserves a gradient structure in the initial refer- 
ence frame (if (j = 0) or in a rotating frame (if s = 0). The parameters uj and s have 
very different roles. One can think, heuristically, that s determines the long time shape 
of the solution while uj affects the rotational behavior of the solution. 

In general though, when both to and s are non-zero there is no obvious gradient 
structure but some of the gradient dynamics features are still present. For flows with 
arbitrary lo ^ Q and appropriately small s, at high concentration 6, we prove the ex- 
istence of tumbling solutions, that is solutions which are periodic in time in a moving 
frame. For arbitrary a;, small s and small 6, the solutions evolve to a steady state in the 
long time limit. 

Intriguing and challenging questions remain to be addressed in the framework devel- 
oped by this paper. For instance, our analysis offers information about the cases when 
s is small and b is either small or large. We can not offer any information about what 
happens, for small s, at intermediary values of h. This would the range of h that, in the 
absence of the flow, corresponds to the transition from isotropic to nematics, from one 
steady state to two steady states. The case of large s seems to be significantly more 
difficult. 

It is of interest to determine how general is the behavior described above. Indeed, 
for quite general ID local nonlinear parabolic equations with bounded trajectories it 
is known that the solutions evolve in the long time limit to either a steady state or a 
periodic solution ([ini,[2Zl). The presence of a non-local term may allow for significantly 
more complicated dynamics (|llj). We do not know if this is the case in here or one can 
prove a Poincare-Bendixson theorem as in f jlUj^. 

The paper is organized as follows: in the next section we consider a general nonlinear 
Fokker-Planck equation and prove its dissipativity. We recognize that both levels of 
complexity can be put into this general form and this gives us the dissipativity at both 
levels. Moreover, when the nonlinear Fokker-Planck equation has a gradient structure 
we analyze its w-limit sets. We obtain thus information about the case when the flow is 
neglected or the flow is present but irrotational. 

In the last section we consider the case when the flow is present and is externally 
imposed, of constant gradient. If s and b are small we prove the evolution to a steady 
state. The proof also shows the existence of a unique steady state by using arguments 
from the dynamical systems theory. 

Finally, we prove our main theorem: the existence of time periodic solutions in the 
appropriate moving frame for small s and large 5. The proof involves the use of a 
non-standard form of the implicit function theorem in Banach spaces. 

2 The dissipativity 
2.1 The general case 

Consider the nonlinear Fokker-Planck equation: 

dtf + deifV + fdelCf) = deef (5) 



3 



where / depends only on time and 9 G [0, 2tt] with V a potential which depends on 

G [0, 2tt\ and may depend on time i.e. V = V{t, 9). We assume that /(t, •) is periodic 
f{t,0) = /(t, 27r),Vt > and since / represents a probability distribution function we 
assume that for all t > it is positive and of mean one. We also assume that V is 
smooth, periodic in the 9 variable and ICf is defined as in (0). 

We prove that any solution of the equation ^ which starts from a nonnegative, 
mean one, smooth initial data will eventually enter a fixed ball in . Let us remark 
that for the case when V = 0, ICf is the Maier-Saupe potential and even initial data 
the existence of a global attractor in any Sobolev norm, was proved in [H]. 

Theorem 1 (The dissipativity) Assume that there exist M,N such that 

\\deV{t,9)\\Loo(^s^) < M,yt > \deek\ < N, V0 

Let f be a solution of |3p starting from a smooth nonnegative initial data, with 
11/(0)11^1 = 1. We have 

ll/(i)lli.<ll/(0)|li.e-*/^ + C + i (6) 

TT 

\\def{t)\\h < (l|9./(0)||i. + ||/(0)||i. • 6(M + bN)) e-*/^ + (M + bN) (sC + ^) (7) 
with 

where C is a constant that appears in the Gagliardo-Nirenberg inequality 

Proof. The existence of solutions is obtained by standard arguments (see also 
The positivity of the initial data is preserved, by the maximum principle, and since 

1 Jsi f{t,e)d9 = we have that ||/(t,-)llLi =l,Vi >0. 

Let us recall the classical Gagliardo-Nirenberg inequalities ( |2H1 ) which state that for 
I < q,r < oo and j, m integers satisfying < j < m and for / an appropriately smooth, 
periodic, mean zero function 

< (8) 

with 

1 7 ,1 m , 1 — a 
pa r a q 

where d is the spatial dimension. 

Taking d = l,j = 0,p = oo, m = 1, r = 2, q = 1 we obtain 

ll/-/l|L^<c||a,/||^//ll/-/l^/=' 

where we denoted by / the average of /. Using the fact that the norm of / is 1 and 
/ = 27 the last inequality implies 

LOO <T^ + C2V3(i + 119^/1^2) (9) 



4 



Also, using Poincare's inequality and / = ^ we have 

L^<\\f\\L^ + \\f-f\\L^<^ + \\def\\L^ 

y 111 



and then 

||/|li2<^ + 2||5,/||i. (10) 
Multiplying by /, integrating over and by parts we obtain 

\irfl f" + lf f"9oV + ]- [ {deelCf)f = - [ {def? (11) 
2. at jgi I J si I J si Js^ 



Using the hypothesis, the bound (jUJ and ||/||li = 1 we get 



^1 / deVfde\<hdeV\\L^\\f\\L^\\f\\ 
Z J si 2 

< f + C2V3(1 + ||*/||,= )) < f + C2./3) + , 1 

Similarly 



Using the last two bounds in lfTT|) we obtain 



Id 2 ^ M + bN 1 MC 2 >iVg.2 1||-..||2 ...x 

2^II-/IIl2(51) < ^ W+ ^ + ^22?^^ -2\\^et\\L2 (12) 

' ~ ~ > ' 

C/4 

so then using (|lUj) . multiplying by 26*/^, integrating on [0,t] and then multiplying by 
e-*/2 obtain 

On the other hand (|12jl can be rewritten as 

^|ii/iii. + ^ii9./iii.<^-^ii5,/iii. 

Using (jlUj) on the right hand side, multiplying by 2e*/^, integrating on [0,t] and 
multiplying by e~*/^ we get 



Il/(*)lli2 + I r \\def{s)\\l,e(^-'y'ds <2C + - + ||/(0)||i.e-*/^ (13) 
Jo ^ 

which gives us an apriori bound on a time integral involving ||5e/||^2- 

In order to obtain the dissipativity of ||96)/||j;^2 let us differentiate © with respect 
to 9. Denoting dgf = F we obtain an equation for F: 

dtF + de [FV + fdeV + Fde{lCf) + fdeeilCf)] = deeF (14) 
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Multiplying by integrating on 5^ and by parts we have 

I at J si ^ J si Js^ 
+i [ dee{lCf)F^ - I fdeeilCf)deF = - [ {dePf 



Hence 



ld_ 



F\\h= f {deV + dee{lCf))fdeF- \ I {d,v + dee{ICf))F' - [ {deFf 
Js^ 2 J SI J si 

<k\\deV\\L^+b\\deek\\L^) [ f + \ f {deFf 

+k\\deV\\Loo+b\\deek\\Loo) [ F^ - [ (deFf 
2 Js^ Js^ 

^ 2 + 11^ IIl2(51)) 2 

where for the last inequality we denoted P = ||/(0)||^2e~*/^ + C + ^ and used © which 
we have just proved. 

Using the fact that F = fg is mean zero and Poincare's inequality the last inequality 
implies 

^ll^lli^ + < {M + bN){P+\\F\\l,) 

Multiplying by e*/^, integrating on [0,t] and multiplying by e~*/^ we have 

WFiml, < (||F(0)||i.+4(M + 6iV)||/(0)||i.)e-*/^ 
+ (M + bN) (^i{C + + I \\Fis)\\l2e^'~'^/Us^ 

< (lli^(O)lli. + 6(M + 6iV)||/(0)||i.)e-*/' + (M + bN){8C + -) 

TT 

where for the last inequality we used (fTT?|) . □ 

Remark 1 Repeating the procedures described above for higher order derivatives allows 
one to obtain, inductively, that the equation is dissipative in any Sobolev norm. 

Remark 2 One can improve the rate of decay into the absorbing ball at the expense of 
having an absorbing ball of bigger radius. 

Remark 3 Note that equation |3p generates a compact nonlinear semigroup. Indeed, 
consider the mapping S{t) : H C H 5/^1(0, 1) H C H 5'j^i(0, 1) which associates 

to an element f the solution at time t starting from initial data f (we denoted by C the 
cone of nonnegative functions and by S/^i(0, 1) the unit sphere in L^{S^)). Then, as 
= 1,'^t > and ICf is the convolution of a smooth kernel k with f, we have 
that ICf is apriori bounded in any Sobolev norm. Thus the equation can be essentially 
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treated as a semilinear parabolic equation and standard arguments give the existence, 
uniqueness and continuous dependence on the initial data which shows that S{t) is a 
semigroup. The compactness is a consequence of the usual smoothing effect of parabolic 
equations. See for details llSj/ . Ch.3. 



2.2 The gradient case 

In the following we assume that there exists a 27r-periodic function W = W{0) such 
that 

V = deW (15) 

We have then that ^ becomes an equation of gradient type with the energy func- 
tional 

£ = f \ogf.f- \ [ ICf-f-f W-f 

(see also |3]) and 

^ = -f \de{\ogf-ICf-W)\^fde (16) 

We show that the presence of this energy functional is enough for proving that the 
w-limit set of any solution if made of steady states. We show that the cj-limit reduces 
to only one steady state if additional symmetry constraints are imposed. 

We first need some properties of the energy functional: 

Lemma 1 Assume that W G L°°(S^). Then the energy functional £ is bounded from 
below along the solutions and it is locally Lipschitz as a functional from n C into M 
(where C denotes the cone of nonnegative functions). 

Proof. The energy £{f) = f^i / log / — ^ Jgi ICf ■ f — Jgi fW is made of three 
parts: the (negative) Boltzmann entropy J^i / log /, the nonlinear potential contribution 
h /51 ' f ^^'^ linear potential part fgi Wf. 

We have that the nonlinear potential contribution part is bounded in L°° thanks to 
the fact that / > and J^i / = 1. Indeed: 

11/ /c/-/||i^<||/c/|U^||/||ii<Pl|L^(/ ff 

A similar argument works for the linear potential part, using the hypothesis W £ 
L°°(5i). 

On the other hand, the function x log x is bounded from below by — ^ which combined 
with the previous observations gives us the boundedness from below of the energy E{f). 

The fact that the entropy is a locally Lipschitz functional in is a consequence of 
the inequality (see also [TB]^: 

Ixlogx — y logy I < C{\x — yl'^ + \x — y\{x^ + y2)) 
Also, the nonlinear potential part of the energy is locally Lipschitz in norm: 



7 



II J ICf-f-ICg-gh2<\\ jK,f{f-g) + J IC{f - g)g\\L2 

<\\mL^\\f-9\\L^ + mf-9)\\L^\\9\\L2 

< ||/C||ii^ioo ll/llii 11/ - g\\L2 + ||/C||i2^ioo||/ - s-IIlz 115-1^2 
=1 

where we used the fact that k is smooth and thus ||/C||2,i^loo, ||/C||2,2^j;^oo are bounded. 
It is easy to check that the hnear potential part is also locally Lipschitz in L^. 

Thus we have that the entropy and the potential parts of the energy are locally 
Lipschitz in , which finishes the proof of the lemma. □ 

We prove now that the nonlinear Fokker-Planck equation evolves to the steady states 
in the norm. 

Lemma 2 For any nonnegative initial data of the nonlinear Fokker-Planck equation 
(EP satisfying l\15]) we have that the tu-limit set of the corresponding trajectory 

n = f{tn) ^ for some sequence (tn)n.eN C M+} 
contains only steady states. 

Proof. The compactness in H^nC of the semigroup generated by Q ( see Remark^ 
together with the fact that all trajectories decay exponentially into a fixed ball suffice 
for having a connected global attractor in n C (see for instance |17j . p. 39, Thm. 
3.4.6; note that there one has a semigroup defined on a Banach space. The fact that 
the semigroup is invariant with respect to a cone, as we have here, will not affect the 
validity of the quoted result, as one can easily check). 

Take an arbitrary nonnegative initial data and consider the omega limit set associated 
to the trajectory starting from this initial data, 57. Observe that all the elements in Q 
have the same energy. Indeed, the energy is decreasing along the trajectories and it is 
bounded from below which means that there exists a c € M so that lim^^^oo E{f{t)) = c. 
We claim that this implies that the elements of O are the steady states. 

Recall that the w-limit set is an invariant set (see |17j . p. 36) so any trajectory starting 
from an initial data in $7 will stay in $7 and have the same energy c. For an initial data 
in O equation (|16j) implies log / — KLf — W = const, Vt > 0. Indeed, if the right hand 
side of H16|) is negative at some time to then it will be negative on an interval around to 
and this would imply that £{t) < <S(0),Vt > to, which is a contradiction. From 
and logf - JCf -V = const, Vt > we have that ft{t,e) = 0,Vt > 0,6* G [0,27r] i.e. 

f{t,e) = f{o,e),yt>o,e€[o,27T]. □ 

In general we do not know if 17 reduces to only one steady state. However, when 
we have a certain type of additional symmetry constraint then there are only finitely 
many steady states which can exist in 17 and only one of them will actually be in 17 for 
a given trajectory. We show this for the simple case when k = ^ cos(20) so /C = /Cms 
and V = hence © reduces to 

dtf + de{fde}CMsf) = deef (17) 

This equation has been analyzed in [Sj where it was proved that if one starts with 
an even initial data then the evenness of the initial data will be preserved by the flow. 



8 



Therefore one can restrict oneself to studying solutions which have this symmetry, i.e. 
solutions of the form 

f{t,e) = ^ + -E^^,yk{t)cos{2ke) 
zvr vr 

where yk{t) = f^"" cos{2ke)f{t,9)de. 

The normalization of the initial data implies yo = 1 and \yk\ < 1 . Then the nonlinear 
interaction potential becomes ICMs{d,t) = cos(26'). 

In this setting the equation (|17)) can be rewritten in terms of Fourier coefficients as 
an infinite system of ODE's: 

yo = 1 

y'l. + Ak'^yk = bkyiiyk^i -yk+i),k = 1,2,... 
For k = 1 we have 



y'l = yi(-4 + 6yi - byiy2) 

which implies that if the yi{0) ^ then yi will be non-zero for all times. 

It is known that there are only three even around zero possible steady states for the 
equation H17() (see [HI), an isotropic one fi{9) = ^ and two nematic ones fr{b)~{(^) = 

-r(6)cos(2S) +r(i,)coa(29) i T. ■ . 

jj^^e-'-Wcos(29)rfe and fr{b)+{0) = ^'i. ^+rm.o.(2e)de ^^^^^ ^v^) depends on b. It is easy to 
check that yi{fi) = and yi{fr(h)±) 7^ 0- The eveness of initial data, which is propagated 
by the flow, leads thus to only three possible elements in the w-limit set Q.. But i7 must 
be a connected set (see |17j ) thus it will necessarily consist of only one element. 
Summarizing the last observations, we have proved: 

Lemma 3 Consider the equation \n\) . 

//yi(0) 7^ then fit, 9) tends, as t — > oo, to either fr(b)- or fr{b)+- Otherwise, if 
yi(0) = 0, we have f{t,9) fi. 



3 Homogeneous gradient flows 

In this section we assume that the flow is externally imposed and of homogeneous 
gradient, that is: 

u{xi,X2) = {T,j=iu}xi,i:j^iUiXi) 

with {tij}ij=i,2 arbitrary constants. 

We consider the moving frame transformation 

f{t,xi,X2,9) = f{t,xi +t{J:^^^u}xi),X2 + t{i:^^^UiXi),9) 
and then the equation for /, becomes an equation for / 

dtf + de[{io + s cos(2e + a))f + dglCf • /] = deef 

where we used the fact that V as given in can be written V = to + s cos (20 + a) 
with uj = ^ the vorticity of the flow, s > uniquely determined and a uniquely 
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determined modulo 27r. The constants uj, s, a can be expressed in terms of u*-, i,j = 1, 2. 
One can easily check that f{t,6 — a/2) satisfies the equation 

dtf + de[{uj + s cos(20))/ + deJCf ■ f] = deef (18) 

We continue working with the last equation. We start by observing that this is an 
equation of the form ^ and the assumptions of Theorem^aie fulfilled. Thus equation 
(|18|) is dissipative. 

Moreover, we can prove that for arbitrary u, with small enough s and b the solutions 
evolve, in the long time limit, to a steady state. The strategy of the proof also shows 
that in the parameter regime given by assumptions ()19() . (|2U|) below, there exists a unique 
steady state solution of (|T8|) . 

Theorem 2 Assume that u> £M and s,b are small enough so that 

1 > s + bV2^{C + - + e)\\dek\\L^ + hdeekU^ (19) 

TT Z 
l>7s + ^\\d^k\\Lo.+h\\dlk\\Lo. 

+bV2^\\k\\L^ (e + (2s + b\\deek\\Loo){8C + -) + C + -] (20) 

y TT TT J 

for some e > with C defined in Theorem Q 

Then any solution evolves to the unique steady state as t ^ oo. 

Proof. We prove the statement in two steps. First we consider the difference 
between two arbitrary solutions and we show that after a certain time to depending on 
the size of the initial data the two solutions will approach each other at an exponential 
rate. In the second step we use step 1 and a contraction argument to show that in fact 
any solution will have to evolve to a steady state. 

Step 1 Consider the difference between two solutions / and g starting from the initial 
data /(O) respectively 5(0): 

dt{f-g)+u^de{f-g)+de[szos{2e){f-g)+delC{f-gyf+delCgif-g)] = dee{f-g) (21) 
Multiply by / — (7, integrate over and by parts: 

9)-fde{f-g) 
[ deif-gf (22) 

I / delC{f-g)-fde{f-g)\<\\delC{f-g)\\L^\\fh2\\deif-g)h2 
< \\dglC\\L2^L^\\f - g\\L2Ri\\de{f - g)\\L2 < bV2^\\dek\\L^Ri\\de{f - g)\\l. 



\tJ if-af-sf s\n{2e){f-gf{e)de-f a,/c(/- 

2 at J SI Jgl Jgl 

We have the following bound 
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where Ri = C + ^ + e (see Theorem ^for the definition of C) and the above inequahty 
holds after the time to when the solution starting from initial data /(O) enters the ball 
of radius Ri in (this time to depends only on the size of /(O), see Theorem^. For 
the last inequality we used the fact that f — g is mean zero and Poincare's inequality. 
Also 

^1 J ^d0g{ICg){f - gf\ < ^WdeelCWL^^L^ llgh^Wf - aWh < ^\\deek\\L°°\\f - gh^ 

=1 

Using the above bounds in H22|) we obtain, for t > to 



2^11/ - 9\\h < -\mf - 9)\\h +4f- 9\\h + bV2^R4dek\\L^\\deif - g)\\h 

+^\\deek\\L^\\f - gh^ 

and by assumption (|19|) and Poincare's inequality we have that the difference between 
/ and g will decay exponentially after time to- 

In order to evaluate the difference in the norm we take the derivative of the 
equation (|T8|) with respect to 6, for two solutions / and g. We denote dgf = F, dgg = G 
and then we have an equation for F — G: 



dt{F-G)+d0[{F-G){io + scos{29))-2ssm{2e){f -g) 

+{F - G)delCg + FdelC{f - g) 
+(/ - 9)dee{1Cg) + fdeeW " 9)] = dee{F - G) 

Multiplying by F — G, integrating over 5^ and by parts we obtain 



4s / cos{2e){f - g){F -G)-2s [ sin(20)(F - G)^ 

' V ' V ' 

+ 11 dee{lCg){F - Gf - I FIC{F - G)de{F - G) 

^ V ' ^ V ' 

+ f {F-Gfdee{JCg)+ [ (f - g)dl{ICg){F - G) 



Te Tr 



[ fdelC{F-G)de{F-G) = - [ de{F - Gf (23) 



where we used Q on the last line three lines to commute dg and /C. 
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We bound each term Tj, i = 1, . . . , 8 by Cj||F — G|||2 or Ci\\de{F — G)|||2 with some 
appropriate constants q. 

|ri + r3 + r4 + r6| < (s + 2s + ^ -6119,6^:1^00)11^-^1112 (24) 

We bound T2 and Tj in the same manner. We show only how to bound T2. Using 
an integration by parts twice we have 

\As ! cos{29){f - g)d0{f - g)\ = \As [ sm{2e){f - gf\ 

< Mf - 9\\h < As\\F - GWl^ 

where we used Poincare's inequaUty in the last relation. 
Then 

\T2 + Tr\ < (4. + ^-\\d',kU^)\\F - G\\l, (25) 

The terms and Tg are treated in the same manner. We show again just how to 
bound one of them 

\n\ < \\F\\L2\\ICiF-G)\\L^\\deiF-G)\\L2 
< R2\\IC\y^L^\\F - G\\l2 \\de{F - G)\\l2 < R2b\\k\\L^V2^\\de{F - G)\\l2 

where we used Poincare's inequality in the last relation with the constant R2 = e + {2s + 
b\\dg0k\\Lac){8C + -), for some e > and the relation holds for the time t > ti after 
which F = dgf e 3^2(0, R2) (see Theorem^. 
Hence 

\n + Ts\ < {R2 + Ri)b\\kU^V2^\\dg{F-G)\\l2 (26) 

where the inequality holds for t > t2 = max{ti,to}) where t2 is the time after which 
/ e B{0,Ri),F = dgf G B{0,R2) (see Theorem^. 

Using bounds . . ()26() into l\2'S\i together with Poincare's inequality and as- 
sumption ()2U() we obtain that ||F — G|||2 decays exponentially after time t2. 

Step 2 Consider the ball B{0, R) in H'^ for some R > max{i?i, i22}- Then Theorem 
1 and the first step show that there exists a time ^3 such that for all t > t3 we have 

S{t) ■.B{o,R)nCn Sli (o, i) ^ 5(o, ii) n c n 5^1 (o, i) (27) 

\\S{t)h - S{t)g4H^ < a\\fo - goWm (28) 

for some a < 1, where S{t) denotes the nonlinear semigroup generated by the equation. 
We denoted by C the cone of nonnegative functions and by S^^i (0, 1) the sphere of radius 
1, centered at 0, in L^. 

Let X = B{0, R) nC n Sli (0, 1). Then X with the metric induced by the norm 
is a complete metric space. Define 

T:X^X,Tf = S{a)f 
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for some a G Q, a > ts. The previous arguments show that T thus defined is a contrac- 
tion. 

Denote ToTo- ■ -oT = T". As T is a contraction we have that as n ^ oo, r"/o fi 
where /i is a fixed point of T. 

Similarly, taking Uf = 5(5)/ for 6 G M - Q,5 > tg, we have U : X ^ X is 
a contraction. Reasoning as before we obtain the existence of a /2 G X such that 

Uf2 = h- 

From Step 1 we have 

lim \\S{na)fi - 5'(na)/2||j^i = 

n— too 

But S{na)fi = T"/i = /i so the last limit becomes: 

ll/i - 5'(na)/2||Hi ^ 0, as n ^ oo (29) 

Recall Hurwitz's theorem in number theory (see for instance [30]) which states that 
for 7 G M — Q there are infinitely many rationals | such that 

An easy consequence is that there exist two sequences {mk)k£N, ink)keN]^k,nk > 
/c, V/c G N such that 

InfcO - mkb\ < i 

Let e > 0. As S{t)f2 is continuous at t = there exists a 6 > such that 

\t\<6^\\Sit)f2-f2\\m<^ (30) 
Thus, for k large enough so that \nj^a — mfc6| < ^ < 5 we have 



\\S{nka)f2 - S{mkh)f2\\m = \\S{nka - mkb)S{mkb)f2 - S{mkh)f2\\m 

= \\S{nka- mkb)f2- f2\\m < - (31) 

where for the first equality we used the semigroup property and for the last inequality 
relation (PU)) . 

On the other hand from ()29() we know that for there exists a rank uq such that for 
Uk > no 

||/i-5(n,a)/2|bi<| (32) 
Putting together and (|H2j) we obtain 

ll/i -/2II//1 < e 

and since e is arbitrary we have that /i = /2- 

We show now that S{t)fi is periodic of arbitrarily small period. 
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Take mk,nk such that {n^a — mkb\ < ^. Assuming without loss of generahty that 
n^a < mkb we have 

S{mkh - nka)fi = S{mkb - nka)S{nka)fi = S{mkb)fi = S{mkb)f2 = f2 = h 

and thus S{t)fi has time period < nikb — n^a < ^. It is well known that a continuous 
function of arbitrarily small periods must be constant. Thus /i is a steady state. 

As the difference of any two solutions tends to as t — > oo this shows that all 
solutions tend to the steady state /i, which must be unique. □ 

Following the proof, one can easily see that we also have 

Corollary 1 // Uy\) and \2U\) are satisfied there exists a unique steady state solution of 
the equation M^) . for an arbitrary smooth potential k satisfying 

Let us observe now that the presence of the flow introduces in the equation a term of 
the form dg[{uj + s cos{26)f)]. The parameters u and s play very different roles. In the 
case when either w or s is zero we have that the equation still has a gradient structure, 
in an appropriate reference frame. 

Lemma 4 Consider equation U8\) . The following holds: 

(i) Ifuj = 0, s 7^ equation hlt^) is of gradient type and the iv-limit set of any solution 
consists of steady states solutions of il^). 

(a) If s = 0,uj ^ make the rotating frame transformation 



Then f satisfies the equation ^ (with V = 0) which is an equation of gradient type. 
Moreover, we have a time periodic solution (in the moving frame) of il^) . namely 
g{9 — ujt) where g is a steady state solution of Ui^) for uj = s = 0. 

(Hi) (an isotropic-nematic patternj Assume that s = 0, a; 7^ 0. Consider solutions 
of il8\) with Maier-Saupe potential, for which the initial data is even around 0, in the 6 
variable. As t ^ 00 we have 



where g = ^ ^/yi(/(0,-)) = and g e {fr{b), f-r{b)} ^/yi(/(0,-)) / (with r{b) a 
function ofb, see J^, Lemma\^and yi{f{0, ■)) = J^^ f{0,xi,X2,0)cos{29)d9). 

Proof, (i) In this case condition (fT5|) is satisfied, as V{6) = cos(20) = ^ sm^26') ^ 
Lemma 121 gives us the conclusion. 

(a) Observe that IC is invariant under the rotating frame transformation 



f{t,e) = f{t,6 + u;t) 






k{e + ujt-e' - iot)f{e' + ujt)de' = {icf){9) 



(33) 



thus / satisfies the equation 



dtf + de{fde{Kf)) = deBf 



(34) 
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which is of gradient type. 

A direct computation, using (1331) , shows that g{9 — ojt) is a solution of (fTH|l. when 
s = (where g \s a, steady state sohition of (|34p '). 

(Hi) Consider the same rotating frame transformation as in the previous part. Using 
the fact that the initial data is the same in the moving frame as in the fixed frame and 
taking into account Lemma |31we are done. □ 

The last lemma suggests that one can think, heuristically, that iv affects the rotational 
behavior of the solution while s affects the long time shape of the solution. The roles 
of LO and s are thus very different. When they are both non-zero one can not reduce 
the equation to a gradient type one just by considering it in the rotating frame. What 
prevents us from repeating the argument above is the fact that in a moving frame 
the linear potential V becomes time dependent. In this situation there is no obvious 
Lyapunov functional. 

Nevertheless, when both s and oo are nonzero, we can prove that for s small enough 
one has time periodic solutions for equation (|18() . Taking into account that (|18|) is 
equation in a moving frame, this implies that in the initial frame we have tumbling 
wave solutions. 

The following argument is done just for the Maier-Saupe potential. This is because 
at the present time only for this potential there is a good understanding of the form of 
the steady state solutions (in the absence of the flow) and of their dependence on the 
concentration intensity parameter b. 

Theorem 3 Consider equation with u; 7^ and let /C = ICms be the Maier-Saupe 
potential. For 6 > 4 and arbitrary < Wi < W2 there is an S depending on Wi and 
W2 such that for s G [0, 5], s < ^{W2 - Wi) and \uj\ € [Wi + s, W2 - s] the equation f7^) 
has a time periodic solution. 

Proof. We present first the strategy of the proof. We are interested in obtaining 
zeroes of a functional ^lA'xMxM— 



where X,y are spaces of functions periodic both in t and 6, whose precise definition 
will not be given because we will see soon that it is more convenient to work with a 
different formulation of the above functional. For that formulation we will make precise 
the functional spaces. 

We will show that for any uj there exists a X^j such that for any A G {—Xuj,Xuj) we 
have 



for some f £ X. This suffices for obtaining the conclusion of the theorem. Indeed, this 
shows that for arbitrary uj there exists an open interval around to, namely (w — A^;, w+Ai^) 
such that for any /i G (u — A^^, a; + A^^), u G {—Xuj, X^j), fJ- — v = lo we have J^{f, /i, z^) = 
for some f £ X. The set {x : Wi < \x\ < W2} is compact so there exists a finite covering 
with intervals of the form (lo — X^,uj + X^j), say 



) = dtf + de [{io + s cos(20))/] + dg [de{JCf) • /] - deef 



J^{f,Lo + X,X) = 



(35) 



{x:Wi< \x\ < W2} = uUi(wfc - X^^,LOk + A,; J 



15 



We take now 5 = mini-^^i iy A^^^. and we obtain the conclusion. 

Returning to (|35j) let us consider the Smoluchowski equation in homogeneous flow 

dtf + de[{LO + A + Acos(2e))/ + de{JCf)f] = deef (36) 
Make the rotating frame transformation f{t, 6) = f{t, 6 + ut). Then H36|) becomes 

dtf + Xdeif + cos(20 + 2u;t)f] + de[de{JCf) • /] = deef (37) 
Let g{6) be an even, nonconstant solution of 

de[dgilCf)f]=dgef 

We know that such a solution exists and it is given by the formula 

gr(6) cos(2e) 
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|^2^gr(6)cos(2e)^0 



for a certain r(b) satisfying a nonlinear equation (see jH]). 

We decompose f{t, 9) = z{t, 0) + g. Then z satisfies the equation 

dtz + Xde[{z + g)+ cos(2e + 2ujt){z + g)] + de[de{lCz)g + de{]Cg)z + de{lCz)z\ = zee 

We prove the existence of time periodic solutions for the above equation, with time 
period ^. Define F : X xR^y 

F{z, A) = dtz + \de[{z + g) + cos(20 + 2iot){z + g)] 

+de[de{JCz)g + dg{ICg)z + dg{ICz)z] - zgg (38) 

with 



X = {z£ ([0,-],H'^[0,2tt]) ,z{t,0) =z(t,2^),Vt G [0, -], 
z(0,e) = z(-,6),dtz(0,9) = dtz(-,9,)\/e £ [0,27r], 

n27r 



/ z(t,^)d0 = 0,Vt G [0,-]} 
Jo ^ 



and 



y = {z£L'[ [0, -], l2[0, 2^]) , z(0, 0) = z(-, 0), V0 G [0, 2^], 



f27r 



z(t,0) = z(t,27r),Vt G [0,-], / z(t,6')d6l = 0,Vt G [0,-]} 

w Jo ^ 

Remark 4 T/ie choice of regularity spaces is somewhat arbitrary, as one can see a 
posteriori that the solution is analytic in time and space. What we need for our proof 
is that the norm of X controls the L°° norm in time and space. We also need that the 
structure of y allows for a simple orthogonal decomposition. 
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We have then that 



F(0,0) = 

We want to apply the imphcit function theorem and obtain the existence of a periodic 
solution for small A. This is a continuation argument, which finds a periodic solution of 
period ^ near one which we already know to exists (for A = 0, see Lemma^. It is due 
to the fact that we are working in a rotating frame that the time periodic solution in 
the initial frame, g{9 — u>t), is stationary in the rotating frame, g in H38|) . 

In order to apply the implicit function theorem we need thus to check that 

Lh = d,F{0, 0)h = dth + deidelCg • h + dgJCh • g] - degh 

as a bounded operator from A' to 3^ is a homeomorphism, i.e., taking into account the 
open mapping theorem, that L is bijective. Nevertheless, this is not the case since, 
as we will see, dim{ker[L)) = dim[range{L)) = 1 so L is a an operator of Fredholm 
index 0. In this situation an implicit function theorem is still possible under a certain 
"non-resonance" condition. This is available for instance in |2j, p. 12. We will present 
it in Lemma below after analyzing the operator L. 

In order to determine the kernel and the range of the operator L we need to study 
equations of the form Lh = f, for / G 3^. Multiplying such an equation by ^^/^g-«'=2ci;t 
and integrating on [0, ^] we have that the equation Lh = f reduces to a decoupled 
system of ordinary differential equations: 

i2ujkhk + dg[de{'ICg)hk + de{K.hk)g\ - dgghk = fk 

where we denote by hk, fk the k-th Fourier mode in time of h, respectively /. 

Thus the problem of determining the kernel and range of L reduces to understanding 
the operator 



L{h) = deeh - dg[dg{ICg)h + dg{lCh)g] (39) 



We have 



Lemma 5 Let L : X ^ Y be a bounded operator as defined in i!^y\) with 

X = {f e h^[o,27t], r f{9)d9 = o}, Y = {f e l^[o,27t], r f{e)de = o}. 

Jo Jo 

Then 

ker(L) = {pdgg,p G M} 
and ^ ^ 

f e Range{L) [ \ [ f{a)da){^ - ^ -^)de = (40) 
Jo Jo 9{G) 27r gif^) 

Thus dim{Ker{L)) = codim{Range{L)) = 1, so L is an operator of Fredholm index 

0. 

Moreover, regarding L as an unbounded operator on Y , with D{L) = X, we have 
that L is a sectorial operator with discrete spectrum contained in the real line. 
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Proof. Standard arguments show that the operator, regarded as an unbounded 
operators on Y , has discrete spectrum and it is sectorial (see for instance ^^jEO])- The 
proof will not be given here. 

Let us observe that the operator can only have real spectrum, i.e. only real eigen- 
values. 

Define 

Ah = --Kh (41) 
9 

and then 



Lh = {he -h^- g{Kh)e)e + (- [de9 - 9{Kg)e\)e = {9{Ah)e)e (42) 

9 9 V ' 

=0 

where the cancellation is due to our choice of g. 

In order to compute the spectrum of L consider the equation 

L{hR + ihi) = {R + iI){hR + ihi) 

with i?, /, /i/j, hj real quantities. 

Separating the real and imaginary parts in the above equation we obtain 

{g{AhR)e)e = RHr - Ihi (43) 



{g{Ahi)e)e = IhR + Rhi (44) 

Multiplying (|43|) by AhR, adding the result to (|44jl multiplied by Ahj, integrating 
over [0, 27r] and by parts we have 

- r g{{AhR)ef - r g{{Ahi)ef = R{AhR, hR) + R{Ahi, hi) (45) 
Jo Jo 



where we used the fact that 



hiAhR= / hRAhi (46) 







Also, let us observe that by multiplying by Ahi, integrating over [0, 27r] and 
by parts, we obtain on the left hand side of the equality the same thing as multiplying 
(|44|1 by AhR^ integrating over [0, 27r] and by parts. This implies the equality of the 
corresponding right hand sides, i.e. 



p2'K r2n r2n f'2n 

/ RhRAhi - / IhiAhi = / IhRAhR + / RhiAhR (47) 
Jo Jo Jo Jo 

and using again we have 

-I{Ahi,hi) = I{hR,AhR) (48) 
which, for 1^0, used in (|45|1 implies 

2lT f2lT 

2 



g{{AhR)er - / g{{Ahi)eY = 
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Taking into account that g > the last equahty imphes that (Ahj)0 = = 0. 

Using this into (|1H|) . (|1H) we obtain hj = = 0. Thus necessarily / = 0, so the 
imaginary part of an eigenvalue must be zero. 

In order to compute the range and the kernel, take f £ Y and denote F{9) = 
Jq f{a)da. Then 

Lh = de{giAh)e) = deF 

implies 

g{Ah)e = F{e) + ci (49) 

In order to determine ci we divide by g on both sides of the last equality and integrate 
on [0, 2tt] obtaining 



= - \ (50) 
Also, integrating on [0, 2tt] both sides of (|49|) we get 

r2Tr 

/ F{e)de + ci • 27r (51) 
Jo 



2-K r'ii^ 
g{Ah)e 



On the other hand 



/•27r pliT 

/ g{Ah)e = - ge-Ah 
Jo Jo 

= 2r(b) [ gi6) s\ni29) [- - -c(h) cos(2e) - -s(h) sm(29)]de 
Jo 9 2 2 

= 2r{bMh)-^-{^-^{l-^j)s{h)]=0 

where we used in the first equality an integration by parts and in the second the relation 
ge = —2r{b)gsm{29){see the definition of g). Also we denoted c{h) = f^^ h{9) cos{29)d9, 
s{h) = J^"" h{9) sm{29)d9 and we used the fact that f^"" g{9) cos{A9)d9 = 1 - f (see 0). 
Using the last computation in 1)511) we obtain 

r2iT 

27rci = - / F{9)d9 (52) 
Jo 

From (|5U() and (|52|) we have that ci is equal to both sides of the equality 

which implies a restriction on F namely ()4U() . 

Observe that relation (|4U|) is the only restriction on the range of L. Indeed, returning 
to (|49() and dividing by g, integrating and recalling the definition of Ah, 1)411) we have 

HO) = ^c{h)g{9) cos(20) + ^s{h)g{9) sm{29) + g{9) f ^^±^da + C2g{9) (53) 
2 2 Jq g{a) 
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where C2 is a constant of integration to be determined. 

Multiplying by cos(2^) and integrating over [0,27r] we obtain 



c{h) = ^c{h)[^ + 1(1 - 1)] + g{e) cos(20) ^^i^j-t^da) dO + C2c{g) 



which implies 

f27r 



c{h){2-^-) = jj g{e)co^{2e){^j^ l^^±^da^de + C2c{g) (54) 

Also, integrating ((SHI) over [0, 27r] and using the fact that we are looking for solutions 
h which have mean zero we get 

-C2 = \c{h)c{g) + g{e) ^^i^±^da^ dO (55) 

Multiplying the last relation by —c{g) and replacing the expression for C2c{g) thus 
obtained into ([5i|l we have 

c{h){2 - \ + \c{gf) = -c{g) j\{e) ^ j' ^^^^^da^ d9 

+ r g{9) cos(2e) ( f l^^l±^d<T\ dB (56) 
Jo \Jo 5(0-) / 

In the last relation we can divide by 2 — | + |c^((?) (which is nonzero as c^(g) > 1 — I, 
for 6 > 4 see |24j . Theorem 2.1) and thus we obtain an expression for c{h) only in terms 
of F and g (as ci can also be determined in terms of F and g, see H50() l. 

Let us check that we can take s{h) to be arbitrary in the representation formula for 
h. Indeed, multiplying ()53() by sin(20) and integrating over [0, 27r] we get 

2/b 

which is always true as 

Jo V ' Jo gicr) ^r{h) Jo g{a) 



-ae/{2r{b)) 



=0 

^w7:^ / J- — dO 

^r{h) Jo g{e) 



=0 



where we used an integration by parts for the equality and (|5Uj) . (|52|) for the cancelations. 

Summarizing: for a given / £ L, satisfying the compatibility condition (|4U() we define 
ci by relation (|5()|) and use this in to obtain an expression for c{h). We use this 
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to get C2 from (|55() and plug everything in ()53() . The h thus obtained wih satisfy the 
equation Lh = f for an arbitrary s{h). 

In particular, if we look for a solution of Lh = we obtain that h £ {pdgg,p E M} 
which gives us the kernel of L. □ 

The properties of L imply that the full operator L has the same kernel as L (where 
now pdgg^p E M is regarded as an element in X). Also, the compatibility condition for 
/ E 3^ to be in the range of L is 

so range{L) has codimension one. 

A technical remark is necessary at this point: we have that for f (z y satisfying the 
appropriate compatibility condition as above there exists some h such that Lh = f. In 
order to make sure that this /i is in A' we need to check the regularity in time. While 
our argument so far does not give but L^ regularity in time, higher regularity in time is 
nevertheless available. This is a consequence of the parabolic nature of the operator L. 

More precisely, the equation 

ik2ujhk + L/ifc = fk 
and the fact that ik2uj is in the resolvent of L, imply that 

hk = {ik2uj + L)~Vfc, V/c E Z - {0} 

with 

\\hk\\H^ < \\{ik20J + L)-^\\L2_,HA\fk\\L^ (58) 

If /c = we have (by the compatibility condition H57|) ') /o E range{L) and thus 

ho = L-^fo 

with 

\\ho\\m < II-^~"^IIl2-.h2||/o||l2 (59) 

As L is sectorial we have that for \k\ large enough \k\ > k^ > the number ik2ui is 
in the resolvent (which we already knew from the lemma) and moreover (see |15j ) 

\\{ik2uj + Z)~^||i2^^2 < C, y\k\ >ko>0 (60) 

where C is a constant independent of k. 

Let C = max{C, \\{ik2u} + L)~^||i2^^2, A; = 0,±1,±2, . . . ,±ko}. We have then 

WhkWn^ < C\\fk\\yk e Z 

which implies that h E L'^{{0, -), H^[0,2Tr]), i.e. the same regularity in time as /. 
Also, using the equation Lh = dth + de[dglCg ■ h + dglCh ■ g] — dggh = / we get dth E 
L2((0, 5),L2(0,27r)) so he C{[0, L^{0,2tt)) (see for instance 0). 

On the other hand we also have that /i is a weak solution of a Cauchy problem for 
the equation Lh = /, with initial data /i(0) E L^(0, 27r). But h{0) = h{^) by time 
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periodicity and then the parabohc regularization effect implies /i(0) = h{^) £ H'^[0, 27r]. 
For / G L2(0, ^,L2(0,27r)) the Cauchy problem with initial data h{0) G H'^[0,27r] has 
a unique solution g £ H^{[0, H'^[0,2tt]) with g{0) = h{0) (see UH]). Moreover it can 
be shown that the uniqueness holds for solutions which are only in C([0, 1/^(0, 27r)). 
Then g = h and thus h will have the necessary regularity in time for being in X. 

The abstract lemma ( which extends the Implicit Function Theorem) that we need 

is 

Lemma 6 fJHJ/, p.l2) Let F : U x V ^ Z with U C X,V C R, where X and Z are 
Banach spaces. Assume that F £ C^{U x V, Z) and: 

• F{xq, Aq) = for some {xq, Aq) £ U x V, Range{DxF{xQ, Aq)) is closed in Z and 

dim{Ker{DxF(xo, Xq))) = codim{Range{DxF(xQ, Xq))) = 1 

• We have the "non-resonance" condition 

DxF{xo, Ao) Range{DxF{xo, Aq)) (61) 
Then there exists a continuously differentiable curve through (2;o,Ao) 

{(x(r), A(r))|r £ {-6, 6), (x(0), A(0)) = (xq, Aq)} 

such that 

F{x{r), A(r)) = 0, for r £ {-5, 5) 

for some 5 > and all the solutions of F{x, X) = in a neighborhood of (xq, Aq) belong 
to the curve specified above. 

The "non-resonance" condition (|6H) in our case is 

- 1^ r g^\a)da)de / 

or 

(2^)V r 9{^)da- r g~\a)da 
Jo Jo 

Recalhng that g{e) = ^''"^'"'j'""' with Z{b) = j^'" e'^^^^^'^'^^^e the last expression 
becomes 

(2vr)2 / {Z{b)f 

Let us observe that this is always true for 6 > 4. Indeed, let h{x) = J^^ e^'^°^^'^^^d9. 
Then 

h'{r{b)) _ 2r 
h{r{b)) ~ T 

(see i). 

Observe that Z{b) = h(r{b)) and multiplying the last equality by ^ we get 
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Z{b) db^ 'b 
and taking into account that r(4) = (see 1^]) we have 

Zib) = 2^e/4 
On the other hand, integrating by parts we have 

=0 

which imphes that Z{b) ^ 27r. 

We want now to check that the solution thus obtained is positive and genuinely time 
dependent (i.e. not a stationary state in the initial reference frame). If A is small enough 
this periodic solution will be near (pointwise in time) 0. More precisely we have 

lkllL°°[0,5]x[0,27r] < C'INIU (62) 

for C independent of n, which is just Morrey's imbedding inequality (see for instance 

ini). 

Thus, for A small enough, z is close to zero in L°° in time and space, and since g 
is positive, z + g will be positive as well. Notice that we have worked in the moving 
frame. Moving back to the non-rotating frame we have that the solution / of the 
equation (|36|) is still time periodic and, pointwise in time, near g{0 + u!t). Using the 
fact that g{9) is nonconstant, this implies that for small enough A this is a genuine time 
dependent periodic solution. Indeed, as g is nonconstant there exists Oi ^ 02 £ [0, vr] 
such that g{9i) / g{02)- Take e = ^\g{92) — ^(^i)!- Then, for A small enough, we have 
\f{t,e) - g{d)\ <e,yee [0,27r],Vt G [0,5] which recalling the definition of / implies 

\f{t,9i+ujt)-g{6,)\<e,yt€[0,-] 

to 

\f{t,e2 + ut)-g{e2)\<e,yte[o,-] 

UJ 

Assume without loss of generality that 62 > Oi,uj > 0. Taking t = < ^ in the 

first relation above and t = in the second one, we obtain a contradiction if we assume 
that / is time independent. □ 
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